library(tidyverse)
library(sf)
library(rmapshaper)
library(here)
library(leaflet)
library(htmltools)
scotland <- read_sf(here("data_raw/SG_NHS_HealthBoards_2019/SG_NHS_HealthBoards_2019.shp"))
class(scotland)
head(scotland, 14)
Simple feature collection with 14 features and 4 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 5512.998 ymin: 530250.8 xmax: 470332 ymax: 1220302
Projected CRS: OSGB 1936 / British National Grid
scotland_geometry <- st_geometry(scotland)

scotland_geometry
scotland_simple <- ms_simplify(scotland)
plot(scotland_simple)

obesity <- read_csv("../data_clean/obesity.csv") %>%
  filter(sex == "All") %>%
  filter(date_code == "2016-2019") %>%
  filter(obesity == "Obese")

scotland_simple %>%
  left_join(obesity, by = c("HBCode" = "feature_code")) %>%
ggplot() +
  geom_sf(aes(group = hb_name, fill = value)) +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), axis.text = element_blank()) +
  ggtitle("Obesity in Scotland") +
  scale_fill_gradient(low = "green", high = "red")
smoking_status <- read_csv("../data_clean/smoking_status.csv") %>%
  filter(sex == "All") %>%
  filter(date_code == "2016-2019") %>%
  filter(smoking_status == "Current smoker")

scotland_simple %>%
  left_join(smoking_status, by = c("HBCode" = "feature_code")) %>%
ggplot() +
  geom_sf(aes(group = hb_name, fill = value)) +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), axis.text = element_blank()) +
  ggtitle("Smoking Status in Scotland") +
  scale_fill_gradient(low = "green", high = "red")
e_cig <- read_csv("../data_clean/e_cig.csv") %>%
  filter(sex == "All") %>%
  filter(date_code == "2016-2019") %>%
  filter(e_cigarette_use == "Currently using")

scotland_simple %>%
  left_join(e_cig, by = c("HBCode" = "feature_code")) %>%
ggplot() +
  geom_sf(aes(group = hb_name, fill = value)) +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), axis.text = element_blank()) +
  ggtitle("E-cigarette use in Scotland") +
  scale_fill_gradient(low = "green", high = "red")
alcohol %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(alcohol_consumption == "Hazardous/Harmful drinker")

alcohol_scotland


  select(value) %>%
  pull()
obesity_scotland <- read_csv("../data_clean/obesity.csv") %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(obesity == "Obese")

obesity_scotland

alcohol %>%
                    filter(sex == "All" | date_code == "2016-2019") %>%
                    filter(hb_name == "Scotland") %>%
                    filter(alcohol_consumption == "Hazardous/Harmful drinker") %>%
                    select(value) %>%
                    pull()
smoking_status_scotland <- read_csv("../data_clean/smoking_status.csv") %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(smoking_status == "Current smoker")

smoking_status_scotland
e_cig_scotland <- read_csv("../data_clean/e_cig.csv") %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(e_cigarette_use == "Currently using")

e_cig_scotland
unique (obesity$sex)

                selectInput("year_selection",
                            "Year Selection",
                            choices = unique(obesity$date_code))
alcohol %>%
                    filter(sex == "All" , date_code == "2016-2019") %>%
                    filter(hb_name == "Scotland") %>%
                    filter(alcohol_consumption == "Hazardous/Harmful drinker")
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkodGlkeXZlcnNlKQpsaWJyYXJ5KHNmKQpsaWJyYXJ5KHJtYXBzaGFwZXIpCmxpYnJhcnkoaGVyZSkKbGlicmFyeShsZWFmbGV0KQpsaWJyYXJ5KGh0bWx0b29scykKYGBgCgpgYGB7cn0Kc2NvdGxhbmQgPC0gcmVhZF9zZihoZXJlKCJkYXRhX3Jhdy9TR19OSFNfSGVhbHRoQm9hcmRzXzIwMTkvU0dfTkhTX0hlYWx0aEJvYXJkc18yMDE5LnNocCIpKQpgYGAKCmBgYHtyfQpjbGFzcyhzY290bGFuZCkKYGBgCgpgYGB7cn0KaGVhZChzY290bGFuZCwgMTQpCmBgYAoKYGBge3J9CnNjb3RsYW5kX2dlb21ldHJ5IDwtIHN0X2dlb21ldHJ5KHNjb3RsYW5kKQoKc2NvdGxhbmRfZ2VvbWV0cnkKYGBgCgpgYGB7cn0Kc2NvdGxhbmRfc2ltcGxlIDwtIG1zX3NpbXBsaWZ5KHNjb3RsYW5kKQpgYGAKCmBgYHtyfQpwbG90KHNjb3RsYW5kX3NpbXBsZSkKYGBgCgoKCmBgYHtyfQphbGNvaG9sIDwtIHJlYWRfY3N2KCIuLi9kYXRhX2NsZWFuL2FsY29ob2wuY3N2IikgJT4lIAogIGZpbHRlcihzZXggPT0gIkFsbCIpICU+JSAKICBmaWx0ZXIoZGF0ZV9jb2RlID09ICIyMDE2LTIwMTkiKSAlPiUgCiAgZmlsdGVyKGFsY29ob2xfY29uc3VtcHRpb24gPT0gIkhhemFyZG91cy9IYXJtZnVsIGRyaW5rZXIiKQoKYWxjb2hvbF9zY290bGFuZF9zaW1wbGUgPC0gc2NvdGxhbmRfc2ltcGxlICU+JSAKICBsZWZ0X2pvaW4oYWxjb2hvbCwgYnkgPSBjKCJIQkNvZGUiID0gImZlYXR1cmVfY29kZSIpKSAKZ2dwbG90KGFsY29ob2xfc2NvdGxhbmRfc2ltcGxlKSArCiAgZ2VvbV9zZihhZXMoZ3JvdXAgPSBoYl9uYW1lLCBmaWxsID0gdmFsdWUpKSArCiAgdGhlbWVfbWluaW1hbCgpICsKICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBheGlzLnRleHQgPSBlbGVtZW50X2JsYW5rKCkpICsKICBnZ3RpdGxlKCJBbGNvaG9sIENvbnN1bXB0aW9uIGluIFNjb3RsYW5kIikgKyAKICBzY2FsZV9maWxsX2dyYWRpZW50KGxvdyA9ICJncmVlbiIsIGhpZ2ggPSAicmVkIikKCmFsY29ob2xfc2NvdGxhbmRfc2ltcGxlICU+JSAKICBsZWFmbGV0KCkgJT4lIAogIGFkZFRpbGVzKCkgJT4lCiAgYWRkTWFya2VycyhsbmcgPSB+bG9uZ2l0dWRlLAogICAgICAgICAgICAgICAgICAgbGF0ID0gfmxhdGl0dWRlLAogICAgICAgICAgICAgcG9wdXAgPSB+SEJOYW1lKQpgYGAKCmBgYHtyfQpvYmVzaXR5IDwtIHJlYWRfY3N2KCIuLi9kYXRhX2NsZWFuL29iZXNpdHkuY3N2IikgJT4lIAogIGZpbHRlcihzZXggPT0gIkFsbCIpICU+JSAKICBmaWx0ZXIoZGF0ZV9jb2RlID09ICIyMDE2LTIwMTkiKSAlPiUgCiAgZmlsdGVyKG9iZXNpdHkgPT0gIk9iZXNlIikKCnNjb3RsYW5kX3NpbXBsZSAlPiUgCiAgbGVmdF9qb2luKG9iZXNpdHksIGJ5ID0gYygiSEJDb2RlIiA9ICJmZWF0dXJlX2NvZGUiKSkgJT4lIApnZ3Bsb3QoKSArCiAgZ2VvbV9zZihhZXMoZ3JvdXAgPSBoYl9uYW1lLCBmaWxsID0gdmFsdWUpKSArCiAgdGhlbWVfbWluaW1hbCgpICsKICB0aGVtZShwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9ibGFuaygpLCBheGlzLnRleHQgPSBlbGVtZW50X2JsYW5rKCkpICsKICBnZ3RpdGxlKCJPYmVzaXR5IGluIFNjb3RsYW5kIikgKyAKICBzY2FsZV9maWxsX2dyYWRpZW50KGxvdyA9ICJncmVlbiIsIGhpZ2ggPSAicmVkIikKYGBgCgpgYGB7cn0Kc21va2luZ19zdGF0dXMgPC0gcmVhZF9jc3YoIi4uL2RhdGFfY2xlYW4vc21va2luZ19zdGF0dXMuY3N2IikgJT4lIAogIGZpbHRlcihzZXggPT0gIkFsbCIpICU+JSAKICBmaWx0ZXIoZGF0ZV9jb2RlID09ICIyMDE2LTIwMTkiKSAlPiUgCiAgZmlsdGVyKHNtb2tpbmdfc3RhdHVzID09ICJDdXJyZW50IHNtb2tlciIpCgpzY290bGFuZF9zaW1wbGUgJT4lIAogIGxlZnRfam9pbihzbW9raW5nX3N0YXR1cywgYnkgPSBjKCJIQkNvZGUiID0gImZlYXR1cmVfY29kZSIpKSAlPiUgCmdncGxvdCgpICsKICBnZW9tX3NmKGFlcyhncm91cCA9IGhiX25hbWUsIGZpbGwgPSB2YWx1ZSkpICsKICB0aGVtZV9taW5pbWFsKCkgKwogIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIGF4aXMudGV4dCA9IGVsZW1lbnRfYmxhbmsoKSkgKwogIGdndGl0bGUoIlNtb2tpbmcgU3RhdHVzIGluIFNjb3RsYW5kIikgKyAKICBzY2FsZV9maWxsX2dyYWRpZW50KGxvdyA9ICJncmVlbiIsIGhpZ2ggPSAicmVkIikKYGBgCgpgYGB7cn0KZV9jaWcgPC0gcmVhZF9jc3YoIi4uL2RhdGFfY2xlYW4vZV9jaWcuY3N2IikgJT4lIAogIGZpbHRlcihzZXggPT0gIkFsbCIpICU+JSAKICBmaWx0ZXIoZGF0ZV9jb2RlID09ICIyMDE2LTIwMTkiKSAlPiUgCiAgZmlsdGVyKGVfY2lnYXJldHRlX3VzZSA9PSAiQ3VycmVudGx5IHVzaW5nIikKCnNjb3RsYW5kX3NpbXBsZSAlPiUgCiAgbGVmdF9qb2luKGVfY2lnLCBieSA9IGMoIkhCQ29kZSIgPSAiZmVhdHVyZV9jb2RlIikpICU+JSAKZ2dwbG90KCkgKwogIGdlb21fc2YoYWVzKGdyb3VwID0gaGJfbmFtZSwgZmlsbCA9IHZhbHVlKSkgKwogIHRoZW1lX21pbmltYWwoKSArCiAgdGhlbWUocGFuZWwuZ3JpZC5tYWpvciA9IGVsZW1lbnRfYmxhbmsoKSwgYXhpcy50ZXh0ID0gZWxlbWVudF9ibGFuaygpKSArCiAgZ2d0aXRsZSgiRS1jaWdhcmV0dGUgdXNlIGluIFNjb3RsYW5kIikgKyAKICBzY2FsZV9maWxsX2dyYWRpZW50KGxvdyA9ICJncmVlbiIsIGhpZ2ggPSAicmVkIikKYGBgCgpgYGB7cn0KYWxjb2hvbCAlPiUgCiAgZmlsdGVyKHNleCA9PSAiQWxsIikgJT4lIAogIGZpbHRlcihoYl9uYW1lID09ICJTY290bGFuZCIpICU+JSAKICAgIGZpbHRlcihkYXRlX2NvZGUgPT0gIjIwMTYtMjAxOSIpICU+JSAKICBmaWx0ZXIoYWxjb2hvbF9jb25zdW1wdGlvbiA9PSAiSGF6YXJkb3VzL0hhcm1mdWwgZHJpbmtlciIpIAoKYWxjb2hvbF9zY290bGFuZAoKCiAgc2VsZWN0KHZhbHVlKSAlPiUgCiAgcHVsbCgpCmBgYAoKYGBge3J9Cm9iZXNpdHlfc2NvdGxhbmQgPC0gcmVhZF9jc3YoIi4uL2RhdGFfY2xlYW4vb2Jlc2l0eS5jc3YiKSAlPiUgCiAgZmlsdGVyKHNleCA9PSAiQWxsIikgJT4lIAogIGZpbHRlcihoYl9uYW1lID09ICJTY290bGFuZCIpICU+JSAKICAgIGZpbHRlcihkYXRlX2NvZGUgPT0gIjIwMTYtMjAxOSIpICU+JSAKICBmaWx0ZXIob2Jlc2l0eSA9PSAiT2Jlc2UiKQoKb2Jlc2l0eV9zY290bGFuZApgYGAKCmBgYHtyfQoKYWxjb2hvbCAlPiUKICAgICAgICAgICAgICAgICAgICBmaWx0ZXIoc2V4ID09ICJBbGwiIHwgZGF0ZV9jb2RlID09ICIyMDE2LTIwMTkiKSAlPiUKICAgICAgICAgICAgICAgICAgICBmaWx0ZXIoaGJfbmFtZSA9PSAiU2NvdGxhbmQiKSAlPiUKICAgICAgICAgICAgICAgICAgICBmaWx0ZXIoYWxjb2hvbF9jb25zdW1wdGlvbiA9PSAiSGF6YXJkb3VzL0hhcm1mdWwgZHJpbmtlciIpICU+JQogICAgICAgICAgICAgICAgICAgIHNlbGVjdCh2YWx1ZSkgJT4lCiAgICAgICAgICAgICAgICAgICAgcHVsbCgpCmBgYAoKCmBgYHtyfQpzbW9raW5nX3N0YXR1c19zY290bGFuZCA8LSByZWFkX2NzdigiLi4vZGF0YV9jbGVhbi9zbW9raW5nX3N0YXR1cy5jc3YiKSAlPiUgCiAgZmlsdGVyKHNleCA9PSAiQWxsIikgJT4lIAogIGZpbHRlcihoYl9uYW1lID09ICJTY290bGFuZCIpICU+JSAKICAgIGZpbHRlcihkYXRlX2NvZGUgPT0gIjIwMTYtMjAxOSIpICU+JSAKICBmaWx0ZXIoc21va2luZ19zdGF0dXMgPT0gIkN1cnJlbnQgc21va2VyIikKCnNtb2tpbmdfc3RhdHVzX3Njb3RsYW5kCmBgYAoKYGBge3J9CmVfY2lnX3Njb3RsYW5kIDwtIHJlYWRfY3N2KCIuLi9kYXRhX2NsZWFuL2VfY2lnLmNzdiIpICU+JSAKICBmaWx0ZXIoc2V4ID09ICJBbGwiKSAlPiUgCiAgZmlsdGVyKGhiX25hbWUgPT0gIlNjb3RsYW5kIikgJT4lIAogICAgZmlsdGVyKGRhdGVfY29kZSA9PSAiMjAxNi0yMDE5IikgJT4lIAogIGZpbHRlcihlX2NpZ2FyZXR0ZV91c2UgPT0gIkN1cnJlbnRseSB1c2luZyIpCgplX2NpZ19zY290bGFuZApgYGAKCmBgYHtyfQp1bmlxdWUgKG9iZXNpdHkkc2V4KQogICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICBzZWxlY3RJbnB1dCgieWVhcl9zZWxlY3Rpb24iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIlllYXIgU2VsZWN0aW9uIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNob2ljZXMgPSB1bmlxdWUob2Jlc2l0eSRkYXRlX2NvZGUpKQoKYGBgCgpgYGB7cn0KYWxjb2hvbCAlPiUKICAgICAgICAgICAgICAgICAgICBmaWx0ZXIoc2V4ID09ICJBbGwiIHwgZGF0ZV9jb2RlID09ICIyMDE2LTIwMTkiKSAlPiUKICAgICAgICAgICAgICAgICAgICBmaWx0ZXIoaGJfbmFtZSA9PSAiU2NvdGxhbmQiKSAlPiUKICAgICAgICAgICAgICAgICAgICBmaWx0ZXIoYWxjb2hvbF9jb25zdW1wdGlvbiA9PSAiSGF6YXJkb3VzL0hhcm1mdWwgZHJpbmtlciIpICU+JQogICAgICAgICAgICAgICAgICAgIHNlbGVjdCh2YWx1ZSkgJT4lCiAgICAgICAgICAgICAgICAgICAgcHVsbCgpCmBgYAoK